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Abstract 

We find exact and explicit self-similar solutions of the axisymmet- 
ric MHD equations of a self-gravitating polytropic gas. These solu- 
tions are able to describe a flat (uniform density) subsonic internal 
core, contracting homologously, of a collapsing cloud with aspherical 
pressure, aspherical gravitational potential, magnetic field and rota- 
tion. Two of the solutions describing the collapse of a rotating and 
magnetized cloud present outflow at the pole or at the equator. 



1 Introduction 

1.1 Previous work 

The gravitational contraction of a cloud has been extensively stud- 
ied numerically and analytically in the context of star formation and 
molecular cloud core collapse. Most of the analytical studies assume 
a polytropic equation of state: D/Dt(Pp~i) = 0, a simple, spherical, 
cylindrical or disk-like geometry, and look for exact or approximate 
solutions using various analytical techniques (e.g. self-similarity). 
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1.1.1 Hydrodynamical collapse 

Most of the available work adresses the problem of hydrodynami- 
cal spherical collapse (Periston 1966, 1969a, 1969b, Larson 1969, Shu 
1977, Hunter 1977, Goldreich & Weber 1980, Yahil 1983, Munier & 
Feix 1983, Bouquet et al. 1985, Whitworth &: Summers 1985) using 
self-similar approaches. 

Larson (1969) and Penston (1969) first reduced the fluid equations 
of a self-gravitating isothermal gas into ordinary differential equations 
of the radius, r only. They also derived an exact solution, that presents 
uniform density and homologous velocity in the internal part of the 
cloud (see Fig. [j]). At some point, the velocity (relatively to the sim- 
ilarity profile) reaches the sound speed and a bifurcation occurs, the 
density decreases and the solution tends asymptotically to the singular 
isothermal sphere. This result has been demonstrated analytically by 
Whitworth & Summers (1985) in the isothermal case (7 = 1) and by 
Bouquet et al. (1985) for any value of 7. Whitworth & Summers find 
that the central density is quantised, i.e. there are bands of central 
density for which the solution does not go through the sonic surface. 
They claim that the Larson-Penston solution, that belongs to the first 
allowed band, represents the subsonic core of a collapsing cloud and is 
able to describe better the gravitational collapse than other solutions 
found by Hunter (1977) and Shu (1977) that belong to the other bands 
and that are much more centrally peaked. 

These studies leave aside important aspects such as rotation or 
magnetic field (Mouschovias 1978). 

1.1.2 MHD collapse 

Some more recent works investigate the collapse of a magnetized gas, 
making various approximations. Galli & Shu (1993) use a self-similar 
approach and develop a perturbative solution, Nakamura et al. (1999), 
Contopoulos et al. (1998) and Basu (1997) consider the thin disk 
approximation, Chiueh & Chou (1994) and Li (1998a, 1998b) neglect 
the magnetic tension and assume an isotropic magnetic support. A 
notable exception is the work of Low (1992) who was able, using a self- 
similar analysis, to find exact solutions of the collapse of the 7 = 4/3 
magnetized polytrope. Low reduces the MHD equations to the static 
equations modified by a radial static force, that is to say, the Low 
solutions are static in the similarity space. An important restriction 
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of his work is that there is an equilibrium between the non radial 
components of the Lorentz, the gravitational and the pressure forces 
and consequently the orthoradial and azimuthal components of the 
velocity field vanish. Thus, as far as the author knows, and in spite of 
the importance of considering these processes together, there are no 
exact solutions in the literature describing a gravitational collapse of 
a magnetized and rotating cloud. 

1.2 Aim of the paper 

In this paper, we find exact and explicit solutions describing the grav- 
itational collapse or expansion of self-gravitating polytropic gas of 
uniform density including magnetic field, rotation, aspherical thermal 
pressure and aspherical gravitational potential. The main difficulty 
and originality of the present study is that the various fields have a 
non trivial angular dependence. 

The derivation of exact solutions is always of great interest, they 
give explicit examples of how the different terms counterbalance and 
can constitute the starting point of further analytical studies like sta- 
bility or bifurcation (Goldreich & Weber 1980, Bouquet el al. 1985, 
Blottiau et al. 1988, Ringeval & Bouquet 2000). Such solutions can 
also be used as benchmarks, crucial in the testing of the complex MHD 
codes developed for studies related to the gravitational collapse. In- 
deed, the importance of code testing cannot be over emphasized. 

In Sect. 2, we apply the dynamical rescaling method introduced 
by Munier & Feix (1983), to the MHD equation of a self-gravitating 
polytropic gas. Considering axisymmetric geometry, density indepen- 
dent of r and homologous velocity fields, we first reduce the system 
into a set of 8 equations that depend on t, the rescaled time, and 6 
the colatitude only and finally reduce it into 7 ordinary differential 
equations of 6. 

In Sect. 3, we restrict the problem further considering a uniform 
density field and homologous velocity. We find various exact solutions 
of the stationary (with respect to the rescaled time, t) hydrodynamical 
and MHD equations. Some of the solutions with rotation and magnetic 
field have radial velocity fields that become positive (outflow) at the 
pole or at the equator. Sect. 4 presents various solutions, able to 
describe a collapse, with an anisotropic thermal pressure. Depending 
on 7, the polytropic index, some of the solutions present outflows at 
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the pole or at the equator. Sect. 5 concludes the paper. 



2 Reduction of the system 
2.1 The equations 

We consider the ideal MHD equations, in spherical coordinates, of a 
polytropic self- gravitating gas. We thus assume that the magnetic 
field is perfectly coupled to the gas, i.e. we do not take into account 
ambipolar diffusion. In the usual notation, we have: 



d t p + div(pv) = o (i) 

p(d t V + V.gradV) = 

-gradP + /9grad$+ — rotB AB (2) 

Mo 

3 t (Pp~ 7 ) + V.grad(P / 9- 7 ) = (3) 

A$ = -4^G/>, (4) 

divB = 0, (5) 

d t B + rot(V A B) = (6) 

2.2 Magnetic field 

In the following, we consider axisymmetric geometries (d^ = 0) and 
use the divergence free form for the magnetic field (Low 1992): 

B r = -^—d e H 1 (t,r,0), (7) 
r l sin 9 
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B = —-d r H 1 (t,r,9), (8) 
r sm 



B 4 , = ^H 2 (t,r,9). (9) 

where H 2 has been introduced in order to have the same physical 
dimension as H\. 

Eq. [B] is automatically satisfied and Eq. || reduces to the two equa- 
tions: 



dtHi + VrdrHi + -V e dgHi = 0, (10) 
r 



d t H 2 + rd r (-L-V^eHx + -V r H 2 
V r sm r 

f y e # 2 - -L-v+drH?) = 0. (11) 



2.3 Dynamical rescaling 



The study of the system of Eqs. |^-^ is carried out using a method 
which generalizes the usual self-similar technique. Following Munier 
& Feix (1983) and Bouquet et al. (1985), we perform the following 
transformation: 

dt = T 2 (t)dt, 
r = a(t) f, 
P(t,r,9) = Tr(t)P(t,r,e), 
p(t,r,9) = d(t)p(t,r,9), (12) 

*(t,r,0) = /(*)*(*, r,0), 
#i(i,r,0) = h^mfcr^), 
H 2 (t,r,9) = h(t)H 2 (t,?,6). 

This procedure generalizes self-similarity since, the transformed 
fields (e.g. p) are allowed to depend both on rescaled space and time 
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whereas the self-similar solutions depend on rescaled space only. The 
time dependence of the new fields allows (at least numerically) the 
study of the convergence toward the self-similar solutions. This proce- 
dure also allows the study of the dynamical stability of the self-similar 
solution in the rescaled space. 

The derivative against time and space are given by: 

d t = -rd? 
t z 1 a 

d r = -df (13) 



a 



With 



V r 



dr 



Ve = r^, (14) 
at 

V$ = rsm6—*, 
at 

we have: 

dr da(t)r . . . _ a - 

v > = r S=^ (i5 » 

d(j> a 

V <f> = ^SUltf— = -3^, 

After some algebra, one finds that the system of Eqs. [I], ||, ||, ^, 
p7o| , |IT| becomes in spherical coordinates (as the study of the angular 
dependence of the fields is our main goal, we write these equations 
explicitly): 



&ff> + 4^ (r 2 V r p) + ^— Q d e (sin 9 % p) = 
i r z \ / r sin t) v / 



dfVr + VrChVr + ^V d d e V r - - * 
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,ct t J a a z a p a z 

h 2 r A 1 1 fl ~ „ /#o\ 1 



^H 2 ch 1^1 + —r—^(hHi x 



a 6 d p por 2 \r r \ r I sin6 



- - 1~ ~ KFfl V, 2 cot6> 

^ + t/^V0 + -v^y* + - - - 

+2r 2 Ve = — 27-4^ + - 

\a t/ era r p a z r 

/iVl 1 n /sin0~ \ 

a D a pnor^smti \ \ r z J 



drK/, + K^-V^ + -VedeVtf, + — ^ H ^ 

n 2 fa f \ ~ /rY 4 1 1 

+ 2r 2 y = . x 

\a r/ a a p po r smO 



sm.6 r V r 
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+ ^,4-^ (sine ^,6) = -4ttG^p, 



1 r h 



d^H 2 + rd^ 



+^d e V e H 2 
r 



r sin a r 



sm ( 



,h 

'h J 



(23) 



We choose the scaling in Eqs. [L2 in order to cancel explicit time 



dependencies in system of Eqs. 1 



. Thus we impose that: 
K x , r 2 (a/a - t/t) = K 2 , 



t 2 (d/d + Za/ 
r 4 d/a = K 3 , -KT 4 /(a 2 d) = K A , 
fr 4 /a 2 = K 5 , h 2 T 4 /a 6 d = K 6 , 
a 2 d/f = K 7 , T 2 (d/dt{7rd-^))/7rd-~t 
T 2 h/h = Kg , 



(24) 



where Ki are real numbers. Although there are 6 variables only and 
nine constraints, the system of Eqs. M admits the following solution: 



r{t) 
d(t) 
fit) 



^/n ujt + 1 , a(t) = (n ujt + 1) A 
(n ujt + l)~ 2 , ir(t) = (Q ujt + 1 



,2A-4 



(25) 



(tt UJt + 1 



i2A-2 



, h(t) = (n UJt + 1 



i3A-2 



where u and A are two real numbers, po is the density at origin and 
the Jeans frequency: f2o = i^Gpo) 1 ^ 2 has been introduced for later 
simplification. 

With Eqs. 24, one obtains the dependence of the nine constants 
Kj on u and A: 




1/2) 



Q w(3A - 2) , K 2 = O w(A 
(n uj) 2 X{X - 1) , K A = 1 , 
1 , K 6 = 1 , K 7 = 1 , 
2O o;(A + 7 - 2) , K 9 = O ^(3A 



(26) 



With Eqs. |24j and |26|, the system of Eqs. p.q-|23| is independent of 
the variable t. It can be seen that due to the transformation stated 
by Eqs. 12 new terms, that can be seen e.g. as source terms, appear 
in the new frame defined by the variable t and r (see e.g. Bouquet et 
al. 1985 for a discussion). 

If lo is positive, the transformation is defined for any value of t 
and describes an expansion, and if u is negative the density becomes 
infinite at t = — 1/(Oqo;) and the transformation describes a collapse. 
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The transformation performed in this paper is slightly more general 
than the transformation performed by previous authors (Munier & 
Feix 1983, Bouquet et al. 1985) who assume that 7r(i) = cf (t), i.e. 
they impose that energy is conserved in the similarity space. This 
choice leads to A = 2 — 7, whereas in the present study A is a free 
parameter, energy is not conserved in the similarity space since a 



source term appears in Eq. 20 



2.4 Reduction into systems of ordinary differ- 
ential equations 

In order to reduce further Eqs. 16-2^, it is worthwhile to reduce these 
equations into systems of ordinary differential equations. 

2.4.1 Ordinary equations of the radius 

The first and obvious possible reduction is to consider fields indepen- 
dent on 6. In this case, Vq, V^, Hi, H2 must vanish and with A = 2 — 7, 
one obtains the following equations: 

Op + ja 9 ? ( f2 ft p) = - fi °^ 4 " ( 2? ) 

dfi + Vrd^Vr + fi ^(3 - 2rf)V r + (O w) 2 (2-7)(l- 7 )f 

d-P 

= ~ ^ + ^$, (28) 
P 



^(f$) = -4vrGp. (29) 

which is the system obtained by Yahil (1983) and by Bouquet et al. 
(1985). In the isothermal case (7 = 1) it reduces to the system ob- 
tained by Penston (1969) and Larson (1969) and extensively studied 
by Shu (1977), Hunter(1977) and Whitworth & Summers (1985). 

2.4.2 Ordinary equations of the colatitude 



Another possibility is to reduce the system of Eqs. 16-23, into a system 
that depends on the angle 9 only. 
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This reduction is obtained by considering the fields: 



t = t/fio , V r (t,r,6) = Sl Q (V(t,6) - Xu)r, 

V e (t, r, 9) = n U(t,0)r , %(t, r, 9) = Q W(t, 9)? . 

P(t,r,0) = p n 2 IL(t,0)r 2 , p(t,r,0) = Po R{t,9) , 

$(i,f,0)=^M)f 2 , 

Hi(t,r,6) = ^Jl^^h^tJ)? 3 , 

H 2 (t,r,0) = ^J^n h 2 (t,9)r 3 . 



(30) 



With the definition of V{9) (first equation), it is easily seen with 
Eqs. H and || that V r (r,9) oc V{9). We also have V e (r,9) oc U(9) 
and V${r,6) oc W{9). 

As for the Larson-Penston solution, these fields diverge when r 
goes to infinity and the solutions are valid in a finite domain only 
(until the velocity reaches the velocity of the sound or magneto-sonic 
waves). 

Another important restriction is due to the fact that the magnetic 
field (see Eqs. [7|, || and ^) is proportional to the radius, r, and van- 
ishes at the origin. Thus, the magnetic pressure compresses the gas 
and enhances the condensation process instead of supporting the gas 
against the gravitational collapse. However, in the following we will 
be able to add a uniform field and to avoid this restriction for some 
of the solutions (see Eqs. 58-59). 



With Eqs. E4L UH and BO, the system of Eqs. 16-23 becomes: 



drR + 3VR + -^—:de (sin OUR) = 2uj R, (31) 
1 snxv 



drV + V 2 + Ud e V - U 2 - W 2 - u V = -2- + 2<j> 
1 R 

"5 ( 2h 2 + -^T-A + - f-V^i)) > (32) 

R \ sm 6 sint) \sm9 // 



d~U + 2VU + Ud g U - W 2 cot 9 - uj U 



R 



+ 
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1 f 1 6 
S \ ~ a (sin 9h 2 ) + —^dgh^x (33) 
R \sm9 snr 9 

+^—^d e hid e (-^—rdehi 
sm \ sm 



d^W + 2VW + Ud e W + UW cot 9 - u W = 

1/3 2 \ 

S —^T7{hide (sin9h 2 ) - -^—p.d e h x h 2 , (34) 

R Vsin 9 smu J 

d^IUT 7 ) + 2VUR-'' + Ud e (ILR-T) = 

-2( 7 - 2)u> ILR -7 , (35) 



+ -^3 e (sin#d e 0) = (36) 
sm ft 



Sy»i + W(9)h x (9) + U(9)d e h 1 (9) = 2uj /h(0), (37) 



&A 2 + 3Vh 2 - 3d e (-^sW)hi + d e (Uh 2 ) = 2uh 2 . (38) 
1 \ sm J 

The parameter A, describing the radius rescaling, does not appear 
any more in these equations because of the special dependence on r 
stated by Eqs. [iO] which indeed is scale invariant. 

This set of equations (Eqs. [11-38) is still quite complex and solving 



it completely is a formidable mathematical problem. We will thus look 



in the next sections for particular solutions of the system of Eqs. 31 
p§| . More precisely, we will look for stationary solutions with respect 
to the t variable. Such solutions are not stationary with respect to 
the physical time t and in spite of clear physical limitations, ought to 
describe an aspherical collapse or expansion with complex processes, 
like magnetic field or rotation included. 
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3 Exact solutions 



In this section, we derive exact solutions of Eqs. ^T|- |3^ stationary 
with respect to t, (i.e. d? = 0) first in the hydrodynamical (i.e. un- 
magnetized) case: hi = hi = and then consider the MHD case. In 
appendix A, the solutions are summarized and given explicitly in the 
physical space. 



3.1 Uniform density field: gravitational po- 
tential and thermal pressure 

Eqs. |3l|, [3^, ^3], |35|, and |36| cannot be easily integrated and in order 
to further reduce these equations, we look for solutions with uniform 
density: 

R(9) = 1. (39) 



3.1.1 Gravitational potential 



One can easily obtain a solution of Eq. |36|: 

MO) = — + 0i (3 cos 2 0- 1). (40) 
6 

If 0i 0, the gravitational potential is aspherical whereas the 
density field is constant. This result can be physically understood in 
the following way. The fields: 

Pa(r,e) = j^(a{a + l)A+ 

(a(a + 1) -6)B (3 cos 2 9- 1)) P" 2 , (41) 



$ a (r,0) = -(A + 5(3 cos 2 0- (42) 

are solutions of Eq. ^ (A and B are two real numbers). The physical 
requirement p a > implies that: 

a{a + l)A + 3B(a{a + 1) -6) > , 

a(a + 1)A - B(a(a + 1) - 6) > 0. ^ ' 

In the limit a — > 2, it is equivalent to A > 0. 
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With these definitions, R and are equal to respectively Po 
and $ a _>2/(^o^ 2 )- I n the limit a — ► 2, the colatitude dependent part 
of /9 Q vanishes whereas this is not the case for This means, that 
in the limit a — > 2, the colatitude dependent part of the density field 
stated by Eqs. ^ and 42 is dynamically negligible but not gravita- 



tionally negligible. 

An explanation for the origin of the aspherical potential will be 
presented in Sect. 3.2.2| . 



If 0i < -1/6 then 0(0) < and 0(7r/2) > 0. The gravitational 
force attracts towards the origin the gas at the pole and moves away 
from the origin the gas at the equator. 

If 0i > 1/12 then 0(0) > and 0(vr/2) < 0. The gravitational 
force attracts towards the origin the gas at the equator and moves 
away from the origin the gas at the pole. 

3.1.2 Thermal pressure 

The thermal pressure given by Eqs. [3(] is proportional to H(9)r, thus 
a negative value of 11(6?) could appear unphysical (negative thermal 
pressure). However with the assumption stated by Eq. [39], it is easily 
seen that: 

P(t, r, 9) = d(tyP + n(t)P(t, f, 9) , (44) 

where Pq is a real number, is also a solution of the problem (<i(i) 7 Po(d(t)R)~' y 
is a constant). Thus choosing Pq > 0, allows to have a thermal pres- 
sure positive until a finite value of r. It still becomes negative for an 
arbitrary high value of r, but with the assumption stated by Eqs. 
as we already said, these solutions are valid only in a finite domain. 
The value of Po can thus be high enough for the thermal pressure to 
be positive in the whole domain of validity. 

With the uniform density field assumption, P / implies that the 
temperature is not a function of the density only as it is usually as- 
sumed. An anisotropic temperature can be due to anisotropic heating 
(Nelson & Langer 1997 who consider an anisotropic UV heating), or 
to isotropic heating and anisotropic optical opacity resulting from the 



anisotropic shape of the cloud (see the discussion of Sect. |3.2,2| on the 
anisotropic sonic surface). 

In Sect. 3 we will consider only the case where P = 0. This implies 
that the thermal pressure does not have any dynamical effect since it 
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is uniform. However, as for the Larson-Penston solution, it is expected 
that at the sonic point, i.e. when the velocity is equal to the sound 
speed, a bifurcation occurs (see Fig. |l|). Thus, thermal pressure plays 
an important role since it induces the bifurcation and the solutions 
are not equivalent to the solutions describing the collapse of a cold 
cosmological cloud that has no thermal pressure (Lin et al. 1965). 
In Sect. 4, we will consider the case where 



3.2 Hydrodynamical equations 




Figure 1: Density and velocity fields corresponding to the Larson-Penston 
solution. The density is constant for r < 3 and the velocity is homologous. 
A bifurcation occurs at the sonic point located at r = 3. 



In this section, we consider the aspherical hydrodynamical case 
when the rotation is non zero and the gravitational potential depends 
on 9 (4>i 0). We look for solutions with uniform density (Eq. ^) 
and assume that: 



V(6) = V -Ui cos 2 9 , U{9) = U x cos 6 sin 9 , 
W{9) = Wism9. 



(45) 



Eqs. ml B2| B3 and 34 lead to: 



E/i = 3Vb - 2u, 

V$ - Wl - ujV = -1/3 - 20i, 
-2lWi + XJ\ + Wl + uUx = 6c? 
Wi = 0ovV = uj/2. 



(46) 
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The first equation derives from the mass conservation whereas the 
second and the third derive from the radial momentum conservation 
(the orthoradial momentum conservation leads also to the third equa- 
tion). The term in the left-hand side of the second and the third 
equation is the centrifugal force. The fourth equation derive from the 
azimuthal momentum conservation. 



With Eqs. 46, we obtain two solutions stated by Eqs. 48 and 51 



The first one describes the collapse of an aspherical and non rotating 
cloud and the second one the collapse of a rotating and aspherical 
cloud. 



3.2.1 Spherical solution 

We first consider the spherical case with U\ = 0, W\ = and 4>\ = 0. 
We then have: 

V = 2/3 uj , u? = 3/2, (47) 

which indeed is the solution found by Bouquet et al. (1985) and by 
Penston (1969a) and Larson (1969) in the isothermal case. 

Bouquet et al. (1985) and Whitworth & Summers (1985), in the 
case 7 = 1, show that at the sonic point, i.e when V = C s (C s being 



the sound speed), the solution stated by Eq. |47|, becomes unstable and 
a bifurcation occurs. After this point, the density decreases and one 
recovers asymptotically the Chandrasekhar solutions (Chandrasekhar 
1967). 

This is illustrated in Fig. [l] that displays the density field corre- 
sponding to the Larson-Penston solution. The bifurcation occurs at 
r = 3. 



3.2.2 Aspherical gravitational potential 

If W\ = (no rotation) and <f>\ ^ (anisotropic gravitational field), 
the solutions of Eqs. ^ are: 

, , i /~3 1+1201 TT — 10 f3 01 w 

< « = ^o = /!^, (48) 

W V 3 ^1+2401 M 
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Due to an aspherical (see Eqs. [35^, [5], 41, 42) gravitational field the 
collapse or the expansion are aspherical. In the limit <f>\ — * 0, one 



recovers the solution stated by Eq. 47. A solution exist only if 0i > 
-1/24. 

It is likely (at least if 0i -C 1) that as for the Larson-Penston 
solution a bifurcation will occur at the points of the sonic surface 
defined by: V 2 + Vq = C 2 S . With this definition, it is easily seen that 
the sonic surface is an ellipsoid defined by the relation: 

^tf + y*) + Q2 °< yo - Ul)2 *> = l, (49) 

where x = r sin# cos 0, y = rsin#sin0 and z = rcosO. The length 
of its axes are: a = C s /(Qq\Vo\) and b = C s /(Q,q\Vq — J7i|), and the 
aspect ratio is given by: 

a = Vo-Uj = 1 + 24^! 

b V 1 + 60i ' 1 ' 

This fact, indeed leads to an explanation for the origin of the aspherical 
gravitational potential since the sonic surface is not spherical and since 
the gas density outside the sonic surface is not uniform. 

If </>i — > —1/24 or (pi — > oo, then V(0) — > oo which means that the 
sonic surface becomes smaller and tends to the origin. 

We will not consider such bifurcation further even if this rather 
difficult (particularly in two dimensions) question requires an accurate 
and careful investigation at some stage. 

If 01 > then a/b > 1 and the cloud or at least the sonic surface is 
oblate. It is prolate in the other case. The velocity field corresponding 
to the three values cj>\ = 1/12 (oblate cloud with aspect ratio equal to 
2), 0i = (spherical solution), 0i = —1/36 (prolate cloud with aspect 
ratio equal to 5/2) are displayed in Fig. ^. 



3.2.3 Rotation 

We now consider a non vanishing azimuthal velocity. 
The parameters of the solution are: 

' U = ±yfe/3 + 80i - 4W? , 

< 01 = 1/12 ; Ux = -u/2 , ( 51 ) 
t V(tt/2) =V = uj/2 , V(0) = u. 
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The gravitational potential is strongly constrainted: <pi = 1/12, and 
0(0) = 0. The sonic surface is oblate with an aspect ratio of 2, the mi- 
nor axis being parallel to the rotation axis. The frequency, to, decreases 
with increasing rotation velocity and the solution becomes stationary 
for W± = l/\/2) meaning that for this value the cloud is rotationaly 
supported. 



3.3 Magnetohydrodynamics 

In this part, we consider the gravo- magnetic condensation and derive 
exact solutions. 

If <%= 0, Eqs. |l| and |7] lead to: 

K R(0)V (9) sin 9 = (52) 

where K is a real number. 

We consider first a problem without rotation and without magnetic 
toroidal field. 



3.3.1 Magnetic field without rotation 

We consider a magnetized, non rotating cloud and we look for solutions 
stated by Eqs. 45. With Eq. 52, we have: 

hi(9) = KUi cos 6» sin 2 (9. (53) 

The corresponding magnetic components are proportional to: 



B r (x -(3 cos 2 9 - 1) r, 
Be oc 3 cos 9 sin 9 r. 



We have: 



6 



h x + 



1 



■dgh 



0. 



(54) 

(55) 
li the 



Lorentz force vanishes (see Eqs. 32-33]). 

Before to continue, we demonstrate that adding a spatialy uniform 
magnetic field is possible. Let us consider the uniform magnetic field: 



sin u \ sin ( 

which means that rotB = 0. Thus, for the field given by Eq. 



B r (t,r,9) 
B (t,r,9) 



b{t)B cos 9 , 
-b(t)B sin (9. 



(56) 
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Figure 3: Velocity field and magnetic field lines corresponding to Eqs. and 
to Eqs. [58| - |59| (with B = 20 and iTC/i = 1). The two time steps, t = and 
t = 3/(4f2 M), are displayed. 

The evolution of the velocity field is self-similar and it is multiplied by 4 between 
the 2 steps. The evolution of the magnetic field is not self-similar because the 
uniform component and the component stated by Eqs. [54] have different temporal 
evolution. 19 



where B is a magnetic intensity. Here again the Lorentz force vanishes. 
With Eqs. || and for the velocity field stated by Eq. 45, one finds that: 

b(t) = (1 + Qouit)- 2 ^^. (57) 



Consequently, Eqs. [l||6] admit the solutions given by Eqs. |48| (and 
by Eqs. 12, 25, ^0, ^5|) with the magnetic field (since the Lorentz force 



associated to the fields stated by Eqs. 54 and pa vanishes): 



B r (t,r,9) = 



B(l + n ojt)- 2V °/ w cos 8 + 

n (l + n uty 2 K(-3Ux cos 2 9 + U x ) r , (58) 



-2V /ui 



sin 9 + 



B e (t,r,9) = -B(l + n ut) 

^//xopo x ^o(l + &oujt)~ 2 K (3Ui cos9s'm9) r 



(59) 



for any value of B and K. In the limit U\ 
Eqs. E8| and 7 — > 4/3 in Eqs. 87 and |88|) and K 



(reached for <pi — > in 
00, one finds that 



and 7 

they also admit the spherical solutions (Eqs. 

Fig. ^displays the velocity field and magnetic field lines stated 
by Eqs. [I?] with the magnetic fields stated by Eqs. [5^-|59| {4>\ = 0, 
B = 20, KUx = 1) at two times steps (t = and t = -3/(4O w)). 
The behaviour of the magnetic field is not self-similar because the two 
components of Eqs. 58-3^ are differentially rescaled. 



3.3.2 Rotation and magnetic field 

In this section we consider a rotating and magnetized cloud. We look 



for solutions stated by Eqs. 45, b3 and: 



h 2 (9) = Q sin 9. (60) 
The corresponding azimuthal magnetic component is proportional to: 

Bs oc sin 9 r. (61) 



Eqs. |3j, |32| |33| and |34J lead to: 



Ux = 3Vb - 2lj, 



V 2 



Wl-uoV Q = -1/3 -20i 



2V Ui + Ul + Wl + uUi 
t VqWi - ujWx/2 = KQU X , 



2Q 2 , 
i + 2Q 2 , 



(62) 
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Figure 4: Stream lines (upper panel) and magnetic field lines (lower panel) of 
the fields corresponding to Eqs. ^3| (B = 0) with parameters corresponding 
to the weak rotation regime (first case). The converging center (r = 0) is 
located in the box centre (x = 15, y = 15, z = 15). The gas is accreted at 
the pole and at the equator. 
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Figure 5: Same as Fig. |] for parameters corresponding to the second case 
(strong rotation): 2QK < W\ < 3QK. The gas is accreted at the equator 
and ejected at the pole. 
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Figure 6: Same as Fig. |] for parameters corresponding to the third case 
(strong rotation): 3QK < W\ < AQK. The gas is accreted at the pole and 
ejected at the equator. 



23 



Figure 7: Same as Fig. |] for parameters corresponding to the fourth case 
(strong rotation): 4QK < W\ < 6QK. The gas is accreted at the equator 
and at the pole. 
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whereas Eq. 38 gives the same constraint as Eq. 31. 

The fourth equation derive from the azimuthal momentum con- 
servation, its right-hand side is the magnetic tension which is propor- 
tional to the product of the poloidal and toroidal field. The term Q 2 
in the right-hand side of the second and the third equations is the 
Lorentz force due to the toroidal field. It tends to compress the gas 
because the magnetic toroidal intensity increases with the distance to 
the z-axis, w. There is no magnetic term due to the poloidal field only 
(proportional to K 2 ) because its contribution vanishes (see Eq. 55). 



As in the non rotating case, one can easily check that adding a uni- 
form magnetic field (Eq. |56|) is still possible (the terms (— sin 9)de(sm OB^) 
cos 9 sin 9d r {rB,p) in Eq. [H] and d r (rV^ cos 6) +dg(V ( f > (— sin 6)) in Eq. || 
cancel). The physical meaning is that a solid body rotation does not 
stretch a uniform field parallel to the rotation axis and an helicoidal 
field (uniform poloidal field along the z-axis and toroidal field propor- 
tional to w) does not induce any azimuthal tension. 

Consequently, one obtains the poloidal magnetic field given by 
Eqs. [58-59 whereas the toroidal field is given by Eq. |6l| . 



The solution of Eq. 62 is given by: 



I (2+8Q 2 -4W 2 )(W 1 -3QK) 2 
10 ~ ±xl (W!-2QK) (Vi -6QK) 



1/ip | : ,2 QKW 1 -2(QKf JT Wi ,,/q (63) 

\r(- /n\ _ V _ Wi/2-2QK T/" ( r\\ _ Wx-2QK 
V(ir/2) -V - Wl ^ 3QK w , V{U) - Wl _ 3QK U- 



Wi, Q and K are free parameters. The transformation W\ — > 
—W\ , QK — > —QK leaves the solution invariant and we can assume 
without restriction that W\ > 0. 



Two cases arise (conditions for to 2 to be positive): 

- Wi < y/1/2 + 2Q' 2 (weak rotation) and W x < 2QK or W\ > 
6QK. The centrifugal force is smaller than the contribution of the 
gravity and the toroidal magnetic forces. 

The gravitational potential is spherical (c^i = 0) if W\+12QK(2Q 2 — 
Wf) = 0. We have V(0)u> > and V(ir/2)u > (inflow). 

If QKUi < 0, the azimuthal magnetic tension tends to decrease the 
azimuthal velocity (see Eqs. ^2|). This is the magnetic braking effect. 
If QKUi > 0, the magnetic tension increases the rotation velocity. 
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This forcing could be related to the Torsional Alfven waves inducing 
collapse studied by Habe et al. (1991). 

This regime is illustrated in Fig. |I] that displays the stream lines 
and magnetic field lines for W\ = 1, Q = 0.52 and K = 1.9. 

- Wi > ^1/2 + 2Q' 2 (fast rotation) and 2QK < W 1 < 6QK. 
The centrifugal force dominates gravity (more precisely the spherical 
part of the gravitational force) and magnetic toroidal pinching. <fii 
is greater than 1/12 and thus gravity tends to move away from the 
origin the gas at the equator (see sect. p. 1.1] ). 

Assuming that, 2QK > \/l/2 + 2Q 2 , three cases have to be dis- 
tinguished: 

- 2QK <W X < 3QK then V(0)w < and V(tt/2)uj > 0, 

uj < (collapse) implies that the gas is ejected at the pole (out- 
flow) and accreted at the equator (inflow). Fig. [| displays the 
stream lines and magnetic field lines for W\ = 2.5, Q = 1.6 and 
K = 0.6. 

As for the weak rotation case with QKU\ < 0, the azimuthal 
magnetic tension tends to brake the clouds. 
Since (f>\ > 1/12, the gravitational potential becomes positive at 
the pole and gravity is responsible of the gas ejection along the 
pole (gravity is the only force along the z-axis). It is not clear if 
this is the same physical mechanism that is responsible of out- 
flows found in numerical simulations of collapse-driven outflow 
(e.g. Tomisaka 1998, 2001). 

- 3QK <W X < AQK then V(0)u > and V{ir/2)u < 0, 

uj < (collapse) implies that the gas is accreted (inflow) at the 
pole and ejected (outflow) at the equator. Fig. || displays the 
stream lines and magnetic field lines for W\ = 3.5, Q = 2.4 and 
K = 0.4, 

Since U\QK > 0, the magnetic tension tends to increase the 
cloud rotation velocity. 

- 4QK < Wi < 6QK then V(0)u > and V(tt/2)uj > 0, u < 

(collapse) implies that the gas is ejected at the pole and at the 
equator (usual behaviour). Fig. ^ displays the stream lines and 
magnetic field lines for W\ = 4.5, Q = 3 and K = 0.3. 
As for the previous case, the magnetic tension tends to increase 
the cloud rotation velocity. 
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In the limit W\ — ► 3QK, one finds the surprising result that uj 
(stationary solution) but, if 18Q 2 K 2 - 4Q 2 - 1 > 0: 



\V(ir/2)\ -> 1/V 6y/18Q 2 K 2 - 4Q 2 - 1 / 0, 

\V(0)\ -> y27 3y/18Q 2 K 2 - 4Q 2 - 1 / 0, (64) 

|C/i| -> v/^VlSQ 2 ^ 2 - 4Q 2 - 1 / 0. 

with V(0)V(7r/2) < and where the sign of V(0) can be positive or 
negative. The physical meaning is that the outflow rate is equal to 
the inflow rate. Indeed the accretion rate, J, is proportional to: 

«r/2 

J(r) = 2 / 2vrr 2 sin 9pV r d6 <xV - Ui/3 (65) 
Jo 



and consequently vanishes (see Eqs. 64). 
Let us summarize the results. 

- We find four different regimes, one corresponds to a low rotation and 

three to a strong rotation (greater than gravity and magnetic 
toroidal forces). For two sets of parameters, one finds outflows 
at the pole or at the equator. 

- In the weak rotation regime with QKU\ < and in the first case 

of the strong rotation regime, the magnetic tension brakes the 
cloud. In the other cases, it tends to accelerate the cloud rota- 
tion. 



Although the magnetic field stated by Eqs. 58, 59 and 61 is not 
force-free, since the azimuthal tension and the magnetic toroidal pinch- 
ing do not vanish, the poloidal magnetic pressure is equal to zero and 
does not support the cloud as it is the case in the numerical sim- 
ulations of gravo-magnetic collapse which generally assume uniform 
initial mass-to-magnetic flux ratio. This fact clearly limits the appli- 
cability of the solutions to collapse for which the magnetic poloidal 
support is not dominant (weak magnetic field or field dominated by 
toroidal component). 



4 Aspherical temperature 

In this section we briefly consider the case of an aspherical thermal 
pressure. We also extend the angular dependence of the various fields 
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and find solutions that present velocity fields not symmetrical with 
respect to the equatorial plan. 

With the uniform density assumption, as we already said in Sect. |3.1.2| , 
this means that the temperature is not a function of the density only 
as it is usually assumed. Although less astrophysically relevant than 
the uniform temperature cases studied in the previous section, such 
temperature gradients are expected when radiative processes are con- 
sidered. 



4.1 Results 

We look for solutions given by: 



' V{9) = V - Ui cos 2 9-U- 
U{ff) = C/i cos sin + C/ 2 sin 

< n(0) = n + iii cos 2 

h\{9) = K(U\ COS Sin V T u 2 sui 

k W(6) = W l sm6 , h 2 (0) = Qsm9. 



T 2 cos sin 
,1a 



~ £i " 1 

1 + II2 cos sin 



(66) 



We then have: 



B r oc (— 3t7i cos 2 + E7i - 3U 2 cos sin 0) r, 
B e oc (3£/i cos sin + 3U 2 sin 2 0) r, 



(67) 



B, 



oc smf r. 



As physically expected, the azimuthal and orthoradial components of 
the magnetic and velocity fields vanish at the pole. 
|33| and 34 give: 



Eqs. 31, 32 



U 



2u, 



V 2 



w? 



9K 2 Ul 



2Q 2 



3^o- 
ujVq 

-2n o -l/3-20a 
-2V Ui + U\ + uTlx + Wl = 

-2n x + 60i + %K 2 XJl + 2Q 2 , 

-2v u 2 + u 2 u x + ujU 2 = -2n 2 - m 2 u x xj 2 , 
n 2 + 6K 2 iw 2 = o, 

VbWi - wWi/2 = KQUx. 

The second, third and fourth equations derive from the radial mo- 
mentum conservation and the third, the fourth and the fifth from the 
orthoradial momentum conservation. The fifth equation shows that 
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there is a non vanishing orthoradial component of the Lorentz force 
at the pole that must be counterbalanced by the thermal pressure. 
Eq. E?5] leads to: 



' (-2C/!)n + ((2 7 - 4)w + 2V - 2£/i)II 1 + 

u 2 n 2 = , 

< (-2U 2 )U + (-2U 2 )U X + (69) 

((27 - 4)w + 2V& - «/i)n 2 = , 
4 ((27 - 4)w + 2y )n + (-c/ 2 )n 2 = 0. 

One obtains the four solutions: 

' n = -iii , (7 - 2)0; + v = , n 2 = 0, u 2 = , 

n = , (7 - 2)uj + v - Ui = , n 2 = 0, u 2 = , 
< u 2 n 2 = 2C/ 1 n , {ul - c/|)n 2 = 2[/ 1 ?7 2 n 1 , (70) 

(7 - 2)u + V = U 1 , 

k n + ni = , t/ 2 n 2 = c^no , 2(7 - 2)w + 2V = u x . 

The first and the second solutions are symmetrical with respect to 
the equatorial plane whereas the third and the fourth are not. 



4.2 Non- symmetrical solutions with respect to 
the equatorial plan 

In this section we give and discuss the two solutions of Eqs. ^38] and 



Eqs. 7C with U 2 7^ and II 2 ^ as they present interesting physical 
features. The solutions symmetrical with respect to the equatorial 
plan are given for completeness in appendix B. 
These two solutions are given by: 



nn + ni 



m 2 ul , u 2 u 2 = 2c/iiio , n 



37-4^2 ) 



l-2iI +4Q 2 



-2Wf 



1/12 - 

1 2-7 
3 4-37 

V{tt/2) = V 



<Pi 
K 2 



2-7 
2( 2 



w T /- 7 ) 2 /8, 
-- 7 w/2 , V(0) = (2 - 



= (3/2 7 



2)u , 



(71) 
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n + u x = o , u 2 u 2 = u ^p , n = -6K 2 ul 

, _ + / -i-n -4Q 2 +2iy 1 2 



y 2(27-3) (6 T -7) 

1 = 1/12 + cj 2 (2 7 - 3)(4 7 - 5)/2, ( 72 ) 
^ = b!^. ^ 1 = 4KQ|E^, C/i = 2(3 7 -4) W , 
V(tt/2) =V = 2(7 - l)w , V(0) = 2(3 - 2j)u. 

For the first solution (Eqs. [7l]), K 2 > => 7 < 4/3, 0i < 1/12 and 
< l/2-n + 2Q 2 . 

The rotation parameter, Wi, is a function of the toroidal magnetic 
intensity, Q, and the polytropic index 7. The rotation tends to prevent 
the condensation whereas the magnetic field, tends to compress the 
gas (see the terms 2W 2 and — 4Q 2 in the expression of to) and it is 
worthwhile to know the global trend of both effects. We have: 

- 4Q 2 + 2W 2 = -2/3 Q 2 37 ^~_ 2 ^" 2 . (73) 

If 7 < (1 + \/7)/3 one has : -4Q 2 + 2W 2 > which means that for 



the solution stated by Eqs. |71|, the global effect of the rotation and the 
magnetic compression is to prevent the collapse. If 7 > (1 + v7)/3 
one has : -4Q 2 + 2Wl < 0. 

For the second solution (Eqs. [72]), one has: K 2 > => 7 < 4/3 or 
7 > 3/2. 
We have: 

- 4Q 2 + 2H^ 2 = -4/3 Q 2 247 ^ 2 ^ 27 (74) 

and 

- 4Q 2 + 2^ 2 > O (13 - V7)/12 < 7 < (13 + V7)/12. (75) 

Thus if (13-^/12 <7< (13+\/7)/12 then -4Q 2 +2VF 1 2 -n -l 
can be positive even with Ho > if Q is large enough, whereas in the 
opposite case, only negative values of Ho (the thermal pressure tends 
to compress the gas) are possible in order for — AQ 2 + 2W 2 — IIo — 1 
to be positive. 
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Two cases arise: 



- —4Q 2 + 2Wf — Hq < 1 (weak magnetic intensity/rotation and/or 
weak aspherical thermal pressure), implies 7/6 < 7 < 4/3. 

One has V(0)w > and V(ir/2)uj > 0. 

- — AQ 2 -\-2Wi — Ilo > 1 (strong magnetic intensity/rotation and/or 
strong aspherical thermal pressure), implies 7 < 7/6 or 7 > 3/2. 

- If 7 < 1, one has V(tt/2)lu < and V(0)u > 0, u < (collapse) 

implies that the gas is accreted (inflow) at the pole and ejected 
(outflow) at the equator. 

- If 1 < 7 < 7/6, to < (collapse) implies that the gas is accreted at 

the pole and at the equator (inflow). 

- If 7 > 3/2, V(tt/2)uj > and V(0)u < 0, u < (collapse) implies 

that the gas is ejected (outflow) at the pole and accreted (inflow) 
at the equator. 



As for Sect. p.3.2| , one finds that the outflow solutions have <f>\ > 



1/12 (positive gravitational potential at the pole). 



5 Conclusion 

In this paper, we apply the transformation proposed by Munier & Feix 
(1983) to the axisymmetric magnetohydrodynamical equations of a 
self-gravitating polytropic gas. We reduce these equations to a system 
of 8 ordinary differential equations of the variable (with Eq. |53| it 
is reduced to 7 equations only). We then derive exact and explicit 
self-similar solutions of these ordinary differential equations capable 
to describe the gravitational collapse of a uniform density cloud with 
aspherical pressure, aspherical gravitational potential, magnetic field 
and rotation. These solutions generalise the Larson-Penston solutions 
in the case where these processes are considered and thus describe the 
subsonic core of a collapsing cloud (Whitworth & Summers 1985). 

For the solutions of the collapse with an aspherical gravitational 
potential, we show that the sonic surface is an ellipsoid that can be 
oblate or prolate. Assuming that a bifurcation occurs at this surface 
(as it is the case for the spherical solution), this suggests that these 
solutions describe the internal part of elliptical clouds. 
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We find that the solutions given by Eqs. Wk pL WJ\ and M de- 



scribing the collapse of a non rotating cloud, are compatible with a 



magnetic field given by Eqs. 58 and [5Sj of any intensity since the mag- 
netic configuration is such that the Lorentz force vanishes. 

The solutions describing the collapse of magnetized and rotating 
clouds with uniform thermal pressure (Eqs. ^) present four regimes, 
one corrresponding to slow rotation and three to rapid rotation. The 
solutions corresponding to these three regimes (rapid rotation) present 
(depending on the rotation and on the magnetic intensity) outflow at 
the pole, outflow at the equator or inflow only. There is a value of 
the rotation (Eqs. |64] ) for which the outflow rate is equal to the inflow 
rate and the solution is stationary (to = 0). 

Four regimes, depending on 7, are also found for one of the (non 
symmetrical with respect to the equatorial plane) solutions describing 
the collapse of a magnetized and rotating cloud with aspherical pres- 
sure (Eq. |7^). One regime (7 < 1) presents outflow at the equator, 
one regime (7 > 3/2) outflow at the pole. 

The two main restrictions of the present work are, first the Lorentz 
force due to the poloidal magnetic field vanishes (if U2 = 0) and thus, 
the magnetic poloidal force is not a support against the gravitational 
collapse (if U2 7^ 0, it compresses the gas), second the solutions diverge 
when r goes to infinity and are valid in a finite domain only. The study 
of the bifurcation that could occur at the sonic (or magneto-sonic) 
point has not been considered here but must be carefully addressed. 
Indeed, a numerical solution of the two dimensional stationary system 
of Eqs. [jJ-H (see e.g. Arthur & Falle 1991 who carry out a similar 
study for a bidimensional self-similar supernovae explosion) must be 
addressed at some stage. It should make possible to avoid the vanish- 
ing poloidal force and should give the solution on the whole domain. 
Such a study should give a sensible hint of the gravo-magnetic con- 
densation. 

Nevertheless, these solutions that generalise the Larson-Penston 
solution have an explicit dependence on the radius r and the colati- 
tude 6 and present a wide class of physical behaviours. They offer an 
opportunity for the testing of numerical codes and a starting point for 
future analytical studies. 
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A Exact and explicit solutions of the 
gravo-magnetic condensation 

In this appendix we summarize the solutions found in this paper and 
give them explicitly. The values of the parameters 7, cj, Uq, Hi, H 2 , 
Ui, U 2 , V , ^i, K and Q being given by Eqs. §@, gjj f§, §]], || ^0 
|S3| , |9l| , p^ , [7l| or [72], po and Pq being free parameters and f2n 
^/4-irGpo, Eqs. admit the solutions given by: 



p(t,r,6) 



Pa 



{l + n Lut) 2, 



(76) 



Ht,r,e) 



o 2 



(1 + fto^) 2 6 



- + 30i cos 2 



>i)r 2 , 



(77) 



P(t,r,9) 



Pa 



+ 



(1 + VLouty-y (1 + n ujt) 4 
(n + ni cos 2 e + n 2 cos 9 sin ( 



(78) 



V r (t,r,9) 

1 + r^o^t 



x (Vq - Ui cos 2 9-U 2 cos 9 sin 0) r, (79) 



Ve(t,r,9) 



On 



1 + On^i 



(C/i cos 6* sin 9 + U 2 sin 2 0) r, 



(80) 



(81) 
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B r (t,r,0) 



(i + n ujt) 2 x 

K (— 3C/i cos 2 O + Ux- 3U 2 cos 9 sin 0) r, 



(82) 



B (t,r,0) 



(i + n ujt) 2 x 

K (3f7i cos 9 sin + 3C/ 2 sin 2 



(83) 



(84) 



With JJi = n 2 = (7,w,IIo,IIi,[/i,Vb, W±,Q and being given 

by Eqs. P. © H H 113' H> S 11' 0' B bein g a free parameter, 
Eqs. [p-gadmit the solutions given by Eqs. ^-81 and 84 and: 



B r (t,r,9) 



B 



[l + fW) 2Vo/ ' 



■ cos 9 + 



(i + n uty 



ET(-3C/i cos 2 9 + Ux) r , (85) 



Be(t,r,0) 



B 



(i + n Lut) 2V o/^ 



sin 9 + 



(1 + n ujt) 2 



K(3U!Cos9sm9) r. 



(86) 



B Aspherical temperature field 

In this appendix we give for completeness the solutions of Eqs. ( 
and |70| that are symmetrical with respect to the equatorial plan (f/ 2 
0, II2 = 0) and briefly discuss their main properties. 
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B.l Aspherical pressure 

Two solutions with W\ = K = Q = are obtained: 

Ho = - IIx , Hp = 3 ^7-5^-1/3+7/4 

- = ±y (T^g^y = (4 - 3 7 )u, , (87) 
F(vr/2) = V = (2 - 7 )£j , V(0) = 2( 7 - 

n = n x = - 1/2 + (37 - 2) (2 - 7 V 2 /4 , 

w = ±^$g*L , ^ = (3 7 _ 4 )o,/2 , (88) 

V(tt/2) = V Q = juj/2 , V(0) = (2 - 7)w. 

These two solutions describe a collapse or an expansion with an as- 
pherical thermal pressure (Fig. displays the velocity field of the so- 



lution stated by Eqs. g§|87| for 7 = 5/3, 7 = 4/3 and 7 = 1.2). The 
gravitational potential can be uniform (0i = 0) or aspherical (0i 7^ 0). 
In the limit <p\ — > and 7 — > 4/3, one recovers the solution stated by 



Eq. 47 



For the first solution (Eqs. pTj ), a solution exists only if 0i < 1/12 
and 1 < 7 < 3/2 or if 0i > 1/12 and 7 < 1 or 7 > 3/2. If 7 < 1 and 
4>\ > 1/12, the solution presents outflows at the pole. 

For the second solution (Eqs. |88|) , one has 4>\ > —1/6 and < 7 < 

2. 



B.2 Aspherical pressure and rotation 

Two solutions with K = Q = are obtained: 

' n = -III , w = ±^4/3 + 801-4^ + 8110 , 
< fa = 1/12 , 7 = 3/2 , Ux = -a//2 , (89) 
, V(tt/2) = Vo = w/2 , V(0) = w. 

n = , cj = ± J4/3 + 8^1 - 41^2 , 
01 = 1/12 + ni/2 , 7 = 1 , E7i = -w/2 , (9°) 
V(7r/2) = V = w/2 , F(0) = w. 

The polytropic index of the first solution (Eqs. |89|) is equal to 3/2 
and the gravitational potential is strongly constrainted (0i = 1/12). 
The solution exists only if n > -1/4 and W\ < y/l/2 + 2U . 
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The polytropic index of the second solution (Eqs. |90|) is equal to 
1. If III = — l/6, the gravitational potential is isotropic {4>\ = 0). The 
solution exists only if W\ < ^/T/W+~2(f)i . 



B.3 Aspherical pressure, rotation and mag- 
netic field 

Two solutions with magnetic fields, rotation, aspherical pressure are 
obtained: 

' n = -ni , 

n = - 1/4-0^+ 1^/2 - ( 7 - 1)(3 7 - 5)^ 2 /2 , 

< » = ±Ji&^,W 1 = 2QK&L, (91) 

U x = (4 - 3 7 )w , 
t V(ir/2) = Vq = (2 - 7 )w , 7(0) = 2( 7 - 

' n = o , 

ni = - l/2-2Q 2 + I^i 2 + (37 - 2)( 7 - 2)w 2 /4 , 

< a, = ± ^ + % iW J + ^ , Wl = QK*=} , (92) 

Ui = (37 - 4)w/2 , 
t 7(vr/2) = 7 = 7^/2 , 7(0) = (2 - 7 )w. 

The first solution (Eqs. |9l|) is defined for 0i < 1/12 and for 1 < 7 < 
3/2 or cpi > 1/12 and 7 < 1 or 7 > 3/2. The second one (Eqs. 
is defined if < 1/3 + 2<pi + 2Q 2 . For both solutions 4>\ can be 
equal to zero (isotropic potential) and the rotation parameter, W%, is 
proportional to QK. 

As for Eqs. |87], if 7 < 1 and <p\ > 1/12, the first solution presents 
an outflow at the pole. 
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